Asymptotic and exact series representations for the 
incomplete Gamma function 



Paolo AmoreJ 

Facultad de Ciencias, Universidad de Colima, 
Bernal Diaz del Castillo 340, Colima, Colima, Mexico 

Abstract. Using a variational approach, two new series representations for the 
incomplete Gamma function are derived: the first is an asymptotic series, which 
contains and improves over the standard asymptotic expansion; the second is a 
uniformly convergent series, completely analytical, which can be used to obtain 
arbitrarily accurate estimates of r(a, x) for any value of a or x. Applications of these 
formulas are discussed. 
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1. Introduction 

This paper is dedicated to the derivation of new series representations for the incomplete 
gamma function, which is given by the integral 



The nice paper by Gautschi pp summarizes many properties of the incomplete 
gamma function and has been used as a valuable reference while writing this paper. 
Since the incomplete gamma function occurrs in a variety of physical and engeneering 
problems, it is very valuable to obtain analytical series representations. Widely used is, 
for example, the asymptotic series j2] 



which can be applied to estimate T for large x. Other series representations, which are 
known in the literature, are valid either for small or large values of x (see The main 
purpose of this paper is to derive new series for the incomplete gamma function, which 
can be used for arbitrary values of x and a; such series could be very useful in dealing 
with a range of applications in physics and engeneering where the incomplete gamma 
function emerges. 

The paper is organized as follows: in section |21 we describe our method and derive 
an asymptotic series for r(a, x), which improves the estimates obtained with eq. ((T]); in 
section El we generalize the results of section |2l obtaining a uniformly convergent series 
representation for r(a,x), which can then be used to estimate r[a,x) with arbitrary 
accuracy. We stress that all the formulas obtained in this paper are fully analytical. 
In section (j^ we consider two applications of the incomplete gamma function, and 
discuss the quality of the approximations obtained. Finally in section (0) we draw our 
conclusions. 

2. The method: asymptotic series 

Techniques such as the Delta Expansion and Variational Perturbation Theory have 
been applied by many authors to physical problems of different nature, ranging from 
classical and quantum mechanics to field theory (see for example jSl 13 El El UD] and 
references therein). In most cases these methods have allowed to obtain very useful and 
precise approximate solutions, which are valid even in the non-perturbative regime, 
where perturbation theory breaks down. The common idea behind these methods 
is simple and appealing: to introduce in the problem an arbitrary parameter, to 
reorganize the expressions by identifying a leading and a perturbative term, both of 
which dependent upon the arbitrary parameter, and finally to apply perturbation theory. 
To a finite perturbative order the dependence upon the parameter survives and its effect 
is minimized by applying the Principle of Minimal Sensitivity (PMS) which allows to 
obtain an optimal value of the parameter. Since the parameter determined in this way 




(1) 





Asymptotic and exact series representations for the incomplete Gamma function 3 



turns out to be a function of the coupling constants in the problem, the method yields 
a truly non-perturbative result, i.e. a result which is not polynomial in the coupling 
constants. 

More recently the author has applied these ideas to obtain a new series 
representation for the Riemann zeta function |1] . The series which was found converges 

oo 

much faster than the standard series = ~t and provides an analytical formula 

n=l 

which is valid for arbitrary s > 0. This formula has also been extended more recently 
to the critical line 

We now come to describe our method: with the purpose of introducing an arbitrary 
parameter we write the term t""'^, appearing in the integral (^, as 

oo r \ 

f^ot^o^a -r)u\{r-uy} ' ^' 



A is the arbitrary parameter. Such expansion is valid as long as 



< 1, i.e. < t < 2A 



A 

and follows directly from the application of the Binomial Theorem. 

Relaxing the constraint of convergence, we use eq. (jH)) inside eq. ((H) and obtain 

r(a, x)^yy — tt(-1)'"" A"^-^-" / e e-* dt . (4) 

Notice that we do not expect eq. (jH) to be an exact series representation of the 
incomplete gamma function, since we have used eq. Q outside the region of convergence. 
For this reason we have used the symbol of approximation in eq. (jH) and we have 
substituted the upper limit in the sum over r, which should be infinite if the series were 
convergent, with a finite integer value R. 

Although the integrals in eq. (0} are still incomplete gamma functions, the exponent 
u is now integer and therefore such integrals can be calculated exactly 



oo 



e e-* dt = e-" y -xP . (5) 
p=o P- 

Finally we can write eq. (0)) as 

- gSS r(.-.)l,t-.): '-^'"' ^' ■ 

Eq. (jni) constitutes a family of asymptotic series for the incomplete Gamma function: 
if such series were convergent the arbitrary parameter A which appears explicitly in the 
equation, would disappear after performing the infinite sum. Although this is not true 
in the present case, given the asymptotic nature of the series, we decide to minimize 
the dependence of the truncated series upon A, as prescribed by the PMS. In for 
example, the PMS allowed to select a series with a faster rate of convergence. 

To first order we obtain the interesting result Xpus = 1 + We thus obtain the 
series 

r(a, x) ^ e-" yyy ^^^^t tt(-I)''"" (l + ■ (7) 
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Figure 1. The ratio Ra{x) as a function of x for fixed values of a. 



Remarkably, to first order, i.e. taking R = 1, eq.© gives the approximate formula 



r(a, x) ^ e ^ (1 + X 



,a-l 



Raix) 



(9) 



In Fig.(P3) our formula (jHI) has been compared with the leading order of eq. (j21) by 
considering the ratio 

Ti{a, x) — r(a, x) 
T2{a, x) — r(a, x) 

where ri{a,x) is given by eq. (jH} and ri{a,x) is given by the leading order term in 
eq. Q. The better behaviour of our formula, eq. (jHI) can be explained noticing that, by 
expanding it for large values of x, we reproduce correctly also the next to leading order 
in eq. (0). 



3. Convergent series 

Although eq. (jH)) provides a useful asymptotic representation for the incomplete gamma 
function, we are interested in obtaining an exact series representation for r{a,x). 
Therefore we need to use the expansion eq. Q only within the region of convergence 
and split the integral in eq. into two pieces: 

- ' - ' (10) 



rfa, x) 



t"-^ e-' dt 



t"-^ e-' dt 



The parameter ^ in this expression, although otherwise arbitrary, must be chosen 
fuUfilling the inequality | < ^ < 2. In the first region one can apply the expansion 
and obtain 



r(a, x) 



r=0 ?i=0 
coo 



+ 



T{a — r)u\{r — u)\ 

oo r 



e e"* dt 



a-l -t 



r=0 M=0 



Via) 



r—u \a—l—u 



T(a — r)u\(r — u)\ 



ul 



+ r(a,eA) 
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As a result, the incomplete gamma function, evaluated at a given point x is obtained 
in terms of the the incomplete gamma function, evaluated at a larger point ^A, which 
lies within a maximum distance of 2A — x. This procedure can now be easily iterated by 
resorting to additional arbitrary parameters, but we prefer to fix A to the optimal value 
found previously, i.e. A = XpMS = 1 + x and to define xq = x, xi = C,X = (,{1 + xq), . . ., 

We then obtain 

r(a, x)=±±± ^ -(-!)— (1 + x,r'-- 

• E^[e-^''a;^-e-«(i+-')(e(l + a;,)r] , (12) 
p=o P- 

a uniformly convergent series representation for r{a,x), which depends upon an 
arbitrary parameter ^. Notice that Xq = ^'^ x + ^^5r- Although one could invoke 
once more the PMS and optimize the series with respect to ^, a large gain in simplicity 
comes from considering ^ = 1, for which one obtains Xq = x + q. 
In this limit we obtain the series 

oo oo r ^ '\ J^(n^ 

g=o r=o ^0 ^0 P! r(a - r) (r - u) ! 

[x + qf -^{l + x + qY^ (13) 

Eq. p3|l constitutes the main result of this paper: when the sums over r and q are 
truncated to a finite order this equation provides an analytical approximation to the 
exact value; the precision of the approximation will be arbitrarily high for sufficiently 
large cutoff values of r and q. We will call the truncated sum r(r^^^,q„^^)(a, x) and define 
the "reduced" function 

p ( \ — (rmax ,qmax) ( T A\ 

where the leading asymptotic behaviour, given by eq. (jH)) has been taken out. 

In the left plot of figure |21 we compare the numerical values obtained for r(7r,x) 
(solid thin line) with the results obtained using the analytical formula corresponding 
to r(6,io)(7r, x). In the right plot we plot the difference between the two curves: the 
maximum error is found for x = 0. 

By taking the limit x — >■ in eq. (fT^ we also obtain the identity 

oo oo r 1 y 1 

1 = y y y _ _ rT(-i)'^-" (i + g)'^-^-"e-^ ( i - - 

,=o.=o«=or(a-r)(r-«)!^ J I ^ ^ 

ooooru-i -1 y-r 

which we have verified numerically§. 

§ Some care should be used because of the condition of convergence < i < 2A. 



p 



Asymptotic and exact series representations for the incomplete Gamma function 6 




Figure 2. Left: Comparison between r{TT,x) and f (6,10) a;); Right: t{Tr,x) — 

f (6, 10) (7^,2:). 



Notice that to improve the accuracy of the series for small x, one could also split 
the integral and write 

T{a,x)= t''-^e-'dt+ t"-^ e'' dt. (15) 

Jx Je 

The first integral could be evaluated with high precision Taylor expanding the 
exponential to a finite order , provided that x < e ~ 1, whereas the second integral 
could be evaluated using once again our formula. 



4. Applications 



A first application of the results obtained in this paper concerns the calculation of the 
probability integral 



P(x) 



e-^y^dt = i erf ^ 
27C Jo 2 \^/2, 



which can be related to the incomplete gamma function by 



erf X 



n 



(16) 



:i7) 



In ^Tj Bagby provides a simple approximation to P{x), given by: 



P(x) — 
^^21 30 



7 e-^"/' + 16 



+ 7 + 



TCX 



1/2 



using the 



In Fig. ini we plot the error over P{x), defined as \P{x) — Papprox\ 
formula of Bagby and our formula eq. (fT^ taking r^ax = 10 and qmax = 20 and relating 
P{x) to r(l/2, x^/2), r(3/2, xV2) and r(5/2, x^/2) respectively. Our analytical formula 
can be made as precise as we wish, both by raising the values of the cutoffs in the series 
and by considering higher orders in the incomplete gamma function. 

Another application of our formula is the calculation of Fresnel integrals 5], which 
naturally occurr in the treatment of diffraction problems: 



C{x) 



77 



COS i —t dt 



S{x) 



sin [ —t dt . 



(19) 
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These integrals can be related to the error function by means of the equation 

C{x) + i S{x) = erf (z) (20) 

where z = ^ (1~0 ^- Our formula already to low order provides an excellent analytical 
approximation to the Fresnel integrals, with somewhat larger errors for small x (see 
Figure IH). Such errors decrease when higher order contributions are considered. The 
dotted lines correspond to the asymptotic formulas: 

+ . (21) 

2 TTX 2 TTX 

5. Conclusions 

In this paper we have considered the incomplete gamma function and obtained two 
new series representations: an asymptotic series, which improves over the standard 
asymptotic expansion of T{a, x), and a convergent series, which to lowest orders can be 
used to obtain arbitrarily precise analytical approximations to the incomplete gamma 
function. As an application of the results obtained in this way, we have used our 
analytical formula to calculate the probability integral and we have compared the results 
with a particular approximation ^1]. By working to order (1, 1) and (5, 5) we also have 
found a very precise analytical approximation to the Fresnel integrals, which occurr 
in the treatment of diffraction. Such approximation could be further improved by 
considering the contributions of higher orders. 

The generality of the ideas upon which the method that we have propesed relies and 
the previous success in dealing with the Riemann and Hurwitz zeta functions motivates 
the effort to obtain similar convergent series representation for other special function. 
Work in this direction is currently in progress. 
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Figure 4. Fresnel integrals C{x) (left plot) and S{x) (right plot) calculated exactly 
(solid line), using the asymptotic formula of eq. (|^ and using eq. (fT!^ to order (1, 1) 
and (5, 5). 
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